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Abstract 

We investigate localization phenomena and stability properties of quasiperiodic os- 
cillations in N degree of freedom Hamiltonian systems and N coupled symplectic 
maps. In particular, we study an example of a parametrically driven Hamiltonian 
lattice with only quartic coupling terms and a system of N coupled standard maps. 
We explore their dynamics using the Generalized Alignment Index ( GALI), which 
constitutes a recently developed numerical method for detecting chaotic orbits in 
many dimensions, estimating the dimensionality of quasiperiodic tori and predicting 
slow diffusion in a way that is faster and more reliable than many other approaches 
known to date. 
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1 Introduction 

In the present paper, we apply the Generalized Alignment Index (GALI) 
method introduced by Skokos et al (2007) [T2] to two types of dynamics, which 
are related to localization phenomena in multi-dimensional Hamiltonian sys- 
tems and symplectic maps. The first type refers to quasiperiodic oscillations in 
the vicinity of what is known as discrete breathers, or exact periodic solutions 
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of multi-dimensional systems, which are exponentially localized in real space. 
The second class of phenomena concerns localization in Fourier space and is 
evidenced by the persistence of quasiperiodic recurrences in the neighborhood 
of normal mode oscillations of nonlinear lattices, observed for example in the 
famous numerical experiments performed by Fermi Pasta and Ulam (FPU) in 
the early 1950's [3ll5lfelfT] . 

In section [2J we present a brief introduction to the GALI method and in 
section [3] we use it to study quasiperiodic motion on tori of low dimensionality 
near a stable discrete breather of a Hamiltonian lattice, with nonlinear on 
site potential and only quartic nearest neighbor interactions. Then in section 
IU we turn to a system of N coupled standard maps and give evidence for 
the existence of both of the above types of localization: First starting with 
initial conditions localized in real space, we find low dimensional quasiperiodic 
motion, which persists for very long times. We also study in this discrete model 
recurrences of its (linear) normal mode oscillations and find that, in contrast 
with the FPU model, the tori associated with them are high-dimensional, 
which may be due to the fact that, unlike the FPU example, each map contains 
on site nonlinear terms, which depend on its individual variables. Finally we 
discuss our conclusions in section [5j 

2 Definition of the GALI method 

Let us consider the 2A r -dimensional phase space of a conservative dynamical 
system, which may be represented by a Hamiltonian flow of N degrees of 
freedom or a 2A r -dimensional system of coupled symplectic maps. In order to 
study whether an orbit is chaotic or not, we examine the asymptotic behavior 
of k initially linearly independent deviations from this orbit, denoted by the 
vectors ~v \,~v i-, k with 2 < k < 2N. Thus, we follow the orbit, using 
Hamilton's equations (or the map equations) of motion and solve in parallel 
the variational equations about this orbit to study the behavior of solutions 
located in its neighborhood. 

The Generalized Alignment Index of order k is a generalization of the Smaller 
Alignment Index (SALI) introduced in [11] and is defined as the norm of the 
wedge (or exterior) product of k associated unit deviation vectors [T2] : 

GALI k (t) =|| h(t) A P 2 {t) A ... A u k {t) || (1) 

representing the volume of the parallelepiped, whose edges are these k vectors. 
We note that the hat (") over a vector denotes that it is of unit magnitude 
and that t represents the continuous or discrete time. 

In the case of a chaotic orbit, all deviation vectors tend to become linearly 
dependent, aligning in the direction of the eigenvector corresponding to the 
maximal Lyapunov exponent and GALL; tends exponentially to zero following 
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the law 

GALI k {t) oc e -l^-^)+(-i—a)+-+(-i--k)]t^ ( 2 ) 
where o\ > ... > are approximations of the first k largest Lyapunov 
exponents of the dynamics. In the case of regular motion, on the other hand, 
all deviation vectors tend to fall on the iV-dimensional tangent space of the 
torus, where the motion is quasiperiodic. Thus, if we start with k < N general 
deviation vectors, these will remain linearly independent on the iV-dimensional 
tangent space of the torus, since there is no particular reason for them to 
become aligned. As a consequence, GALL; in this case remains practically 
constant for k < N . On the other hand, for k > N, GALI^. tends to zero, since 
some deviation vectors will eventually become linearly dependent, following 
power laws that depend on the dimensionality of the torus. According to our 
asymptotic analysis of determinants entering in an expansion of ([1]) in a basis 
of eigenvectors following the motion, one obtains the following formula for the 
GALI fc , associated with quasiperiodic orbits lying on m-dimensional tori [I]: 



GALlJn) oc < 



constant, if 2 < k < in 

j^, ifm<k<2N-m ■ (3) 
^^htj, if2N-m<k<2N 

The computation of such determinants, however, is hardly the most efficient 
way to compute the GALL;, especially in cases where the dimension of the 
system N becomes large. For this reason, we have recently introduced and 
employed a method based on Singular Value Decomposition to show that the 
GALIfc can be very efficiently computed as the product of the singular values 
of a 2N x k matrix [2]. 



3 Applications to the dynamics near discrete breathers 

In this section we examine the dynamics in the vicinity of discrete breathers in 
a one-dimensional Hamiltonian lattice with quartic on site potential and linear 
dispersion terms in its nearest neighbor particle interactions. For this purpose, 
we start with initial conditions near the exact breather and check whether 
the motion remains quasiperiodic or becomes chaotic. We accomplish this by 
computing the GALI indices along the reference orbit. If the breather is stable, 
the GALI method can be used to determine the dimensionality of the tori 
surrounding the breather in the 2iV-dimensional phase space. This dimension 
is generically iV and equals the number of frequencies that are being excited 
in the neighborhood of the exact breather. As for chaotic motion diffusing 
slowly away from these breathers, it is rapidly and efficiently predicted by the 
exponential convergence of all GALI indices to zero. Let us take, for example, 
a Hamiltonian lattice of anharmonic oscillators, which involves only quartic 
coupling terms and hence presents strong localization phenomena due to the 
absence of phonons [7J. More specifically, we study here a system with on site 
potential: V(x) = |[1 — e cos{ojdt)\x 2 — \x A , described by the Hamiltonian [9]: 
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Fig. 1. (a) Initial conditions for the 31 particle lattice, obtained from homoclinic 
intersections of the spatial map <J> n with all velocities set to zero. Here D = 1, 
representing a quasiperiodic breather (Case I in the text), (b) We now perturb the 
initial conditions, so that the amplitudes of the central particle and its 2 neighbors 
are slightly different than those in (a) (Case II in the text). 



Hit) = + ~[1 - ecos{u> d t)\xl - ~4 + j(x n+1 - x n ) 4 }, (4) 

where e, uj d are the amplitude and frequency of the driver respectively. The 
equations of motion of each oscillator clearly are: 

x n = K(x n+ i - x n f + K(x n _i - x n f - [1 - e cos(w d t)]x n + x z n . (5) 

It can be shown that the solutions of equation ([5]) can be written as a product 
of a time-dependent function G(t) and a map $ n , depending only on the 
oscillation site n: 

x n (t) = $ n G(t), (6) 

n — 1, 2, N. Substituting in ([5]), one finds that the equations of motion 
can be exactly separated into a spatial and a time-dependent part each of 
which is equal to an arbitrary constant C. Thus, we arrive at a periodically 
driven Duffing equation: 

G+[l — £ cos(u d t)]G = -CG 3 , (7) 

satisfied by G(t) and a discrete map: 

C$ n + K ($ n+1 - $ n ) 3 + - $n) 3 + € = 0, (8) 

which yields the spatial evolution of the system (jl]). Following [9], we consider 
a lattice of 31 particles, with C = 1, K = l,e = 0.7 and ajj, = 2.6355. The 
initial conditions of the driven Duffing equation are chosen on an exact exact 
periodic orbit with G(0) = D, G(0) = and frequency u;& = Ud/2, if D = 
Df, = ±1.2043. Furthermore, if the initial conditions lie on the homoclinic 
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Fig. 2. Time evolution of GALI/% indices for: (a) D = 1.1 and initial energy 
H(0) = 0.133, see g]), (b) D = 1 and H(0) = 0.097, (c) D = 0.9 and H(0) = 0.07 
and (d) for D = 0.8 and H(0) = 0.048. Notice that in (a) and (b) the breather 
remains quasiperiodic for very long times, lying on a one-dimensional torus, which 
implies that only one main frequency is excited. In (d) the motion is evidently 
chaotic and delocalizes, drifting away from the breather. The slopes of GALI^ evo- 
lution coincide with the laws presented in formulas (|2|) and ([3|). 



intersections of the invariant manifolds of the saddle point at the origin of the 
map, the solution becomes an exact breather with x n (0) = (Fig. [TK), 

which is linearly stable [9] . 

Our first goal, therefore, is to examine the motion near this breather by per- 
turbing the initial conditions in its neighborhood and applying the GALI 
indices to characterize the resulting orbits. In this study, we have made two 
kinds of perturbations: In the first case we change only the factor D, while in 
the second both the factor D and the spatial initial conditions of the map are 
varied. 

Case I: 

Let us start by perturbing the initial conditions G(0) = D (G(0) = 0) away 
from D = Db = 1.2043,, while the initial conditions of $ n are fixed (Fig. [TK). 
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Fig. 3. (a) Oscillations of particles 16, 17, 18 of the lattice show the localization and 
quasiper iodic properties of the motion near the discrete breather. This is verified 
by the time evolution of GALI^ indices in (b) for D = 0.9 and H (0) = 0.09 and in 
(c) for D = 0.8 and H (0) = 0.062, showing the orbits lie on a 3-dimensional torus. 
However, in (d) we take D = 0.7 and Hq = 0.042, perturbing the initial conditions 
as shown in Fig. [TJa. Here, the slopes of the GALI& indices clearly demonstrate that 
the motion chaotically diffuses away. 



In our example, we use a 31 particle Hamiltonian, with C = 1, e = 0.7 and 
uj d = 2.6355. 

Decreasing the value of D < Db, we observe that the oscillations for D — 1.1 
are quasiperiodic and remain so for long integration times, while the particles 
oscillate with frequencies close to that of the exact breather, as shown in Fig. 
|2j This is in agreement with the fact that even GALI2 (representing the area 
of the parallelogram formed by the two unit deviation vectors) tends to zero 
following the power law t~ l . When D becomes 0.9, the GALI^ indices begin to 
fluctuate, indicating that the orbit is near the edge of the corresponding regular 
region in phase space (Fig. [2k). Moreover, near t = 10 6 , GALI 2 begins to tend 
to a constant value, suggesting that the motion has become quasiperiodic with 
3 frequencies. Finally, as D reaches the value 0.8, quasiperiodic recurrences 
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Fig. 4. When we choose initial conditions as shown in Fig. [TJa (Case II), we find 
that with D = 1.2 recurrences break down, while (a) the four largest Lyapunov 
exponents appear to be decreasing to zero. Note in (b) that the difference between 
the two largest Lyapunov exponents, 01 — cr 2 , is nonzero and tends to a constant, 
explaining why in (c) we observe exponential behavior of GALI^,. (d) This agrees 
with the fact that the particles away from the center of the lattice start to oscillate, 
indicating derealization and the breakdown of recurrency. 



break down and the motion becomes chaotic. 
Case II: 

We now perturb slightly the initial conditions of the discrete map $ n (as 
indicated in Fig. db) and observe that the motion now occurs on a 3D torus, 
indicating the appearance of 3 frequencies in the system. This quasiperiodic 
breather behavior lies on a 3D torus, in phase space, as D varies from 0.75 
to 1 (Fig. [3]). In fact, for 1 < D < 1.2, the chaotic character of the orbit is 
revealed by the exponential decay of all GALL; indices. 

For D = 1.2, in fact, recurrences break down and derealization sets in, as 
particles away from the center of the lattice gain energy and begin to oscillate 
(see Fig. HJi). The exponential convergence of GALI indices to zero indicates 
the chaotic behavior of the motion. Fig. shows that there are only three 
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positive Lyapunov exponents, indicating three instability directions about the 
orbit in phase space. The exponents o\ and <j 2 are almost identical, having a 
difference of the order 10 -5 (Figfjb). As a consequence, GALI 2 oc e - ^ 1-0 " 2 - 14 
begins to converge exponentially when t becomes greater than 10 -5 , while the 
other GALIs have reached the value 10 -16 exponentially, long before 10 4 (see 

Fig. IS:). 



4 Application to the dynamics of N coupled standard maps 

The standard map, also referred to as the Chirikov-Taylor map [8], is an area 
preserving map, which appears in many physical problems. It is defined by 
the equations: 

x n +i = x n + y n+1 , (9) 
y n +i = y n + Ksm(x n ), 

where both variables are taken modulo one in the unit square. This map de- 
scribes the motion of a simple mechanical system called a kicked rotator. It 
may be thought of as representing a pendulum rotating on a horizontal fric- 
tionless plane around a fixed axis and being periodically kicked by a nonlinear 
force at the other end, at unit time intervals. The variables x n and y n de- 
termine, respectively, the angular position of the pendulum and its angular 
momentum after the n th kick. The constant K measures the intensity of the 
nonlinear "kicks". 

The standard map describes many other systems occurring in the fields of me- 
chanics, accelerator physics, plasma physics, and solid state physics. However, 
it is also interesting from a fundamental point of view, since it is a very simple 
model of a Hamiltonian system of 2 degrees of freedom that displays order 
chaos. 

For K = 0, the map is linear and only periodic and quasiperiodic orbits are 
allowed (depending on the value of yo), which fill horizontal straight lines in 
the x n ,y n phase plane of the map. However, when K ^ 0, the periodic lines 
break into pairs of isolated points half of which are stable and half unstable 
alternatingly. The stable periodic points possess families of quasiperiodic or- 
bits forming closed curves around them, while the unstable ones are saddles 
with thin chaotic layers close to their invariant manifolds. Moreover, as K > 
grows, the size of the chaotic layers increase and occupy larger and larger 
regions of the available phase plane [SfTO] , see Fig. above. 

Let us consider a system of N coupled standard maps described by the fol- 
lowing equations: 

x n+l = x n + Vn+li (10) 

K 3 

y J n+i=y J n + TT sin(27rx J n ) - — {sin[27r(^ +1 - x 3 n )] + sinpTr^ -1 - x 3 n )}} 

ZlX Z7T 



8 



Fig. 5. Poincare surface of section of a single standard map for K = 2. Note the 
presence of a stable fixed point at (x, y) = (0.5, 0.0) (and its modulo equivalent 
at (x,y) = (0.5,1.0)), surrounded by a large "island" of quasiperiodic orbits. The 
smaller islands embedded in the large chaotic "sea" correspond to a stable periodic 
orbit of period 2. 

with j = 1, N, fixed boundary conditions x = X]y + i = and (3 the coupling 
parameter between neighboring maps. In what follows, we shall examine the 
localization properties of this system comparing the dynamics with what is 
observed when one studies a one-dimensional Hamiltonian system like the 
FPU lattice. 

A. Initial conditions localized in real space 

Let us begin with the question of the existence of stable discrete breathers in 
this model, when the coupling parameter {3 is small. Note that our equations 
have an exact stable fixed point, when all particles are located at (x, y) = 
(0.5,0.0) (see Fig. [5]). Taking N = 20 coupled standard maps, we may, there- 
fore, look for localized oscillations taking as initial condition (Rl): (xj,yj) = 
(0.5,0.0), V j 7^ 11, perturbing only one particle at (xn,yn) = (0.65,0.0) and 
fixing the parameters (3 = 0.001, setting Kj = 2, with j = 1, 20. If a stable 
discrete breather exists, the orbits are expected to oscillate about these initial 
conditions quasiperiodically for very long times. 

Indeed, when we evaluate the GALI& indices for n = 10 6 iterations in Fig. 
[6^ and display their evolution for k = 2, 7 on a logarithmic scale, we find 
that GALI 2 fluctuates around a non-zero value while all the other GALIs de- 
cay to zero following power laws. This implies that the motion is that of a 
quasiperiodic orbit that lies on a 2-dimensional (2D) torus. In fact, for ini- 
tial conditions (R2): (xj,yj) = (0.5,0.0), V j ^ 11,12, perturbing two particles 
Vix) = (0.65,0.0), (xi 2 ,t/i 2 ) = (0.55,0.0) for the same parameters as in 
the Rl experiment, we detected regular motion that lies on a 3D torus! This is 
demonstrated by the evolution of the GALIs in Fig. El where not only GALI 2 
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Fig. 6. (a) GALIfc, k = 2, 7, for for N = 20 - coupled standard maps with Kj = 2, 
(3 = 0.001 and initial conditions Rl (see text). GALI2 fluctuates around a non-zero 
value implying a regular motion that lies on a 2D torus, while the slopes of the 
algebraically decaying indices are in agreement with ([3]).(b) GALI^ for the initial 
conditions R2 (see text) imply that the motion lies on a 3D torus, since GALl2 t 3 
tend to a non-zero number. This quasiperiodic localized dynamics of R2 is clearly 
seen in (c), where we plot the x n oscillations of the 11 th (light gray color), 12 th (gray 
color) and 13 th (black color) map. From the last 10 5 iterations shown in this figure, 
it is evident that the main part of the energy is confined only in the "middle" 3 
maps and is not shared by any of the other degrees of freedom of the system. 



but also GALI3 fluctuates around a non-zero value. 

Such localized regular dynamics becomes evident in Fig. [6)3, where we exhibit 
the oscillations of the x n -coordinate of the 11 th , 12 th and 13 th maps of the 
system for the last 10 5 iterations of the R2 experiment. As seen in the figure, 
the motion is indeed quasiperiodic and confined to the middle 3 maps, as 
all other degrees of freedom do not gain any significant amount of energy. 
Similarly, we have observed that exciting 3 central particles (j = 10, 11, 12) 
gives a quasiperiodic motion on 4D tori. This is reminiscent of the motion near 
discrete breathers of Hamiltonian systems and suggests that the excitation of 
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each particle adds one extra frequency to the motion. 



B. Initial conditions localized in Fourier space 

Finally, let us investigate the phenomenon of FPU recurrences in our system 
of coupled standard maps. To do this, we first need to derive the linear normal 
modes of the model, in order to choose initial conditions that will excite only 
a small number of them. Thus, keeping only the first term in the Taylor 
expansion of the sine-function in ffTOl) . we obtain the following system of 
equations: 



x 
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(12) 



with 7 = 1 + 2/3. 

Using now well-known methods of linear algebra, we diagonalize the matrix 
A writing T> = V~ l AV^ where V is an invertible matrix whose columns are 
the eigenvectors of the matrix A and T> is diagonal, T> = diag[Xi, X 2 n], 
provided the eigenvalues Aj are real and discrete. 

Our case, of course, involves oscillations about a stable equilibrium point and 
hence the above system has 2N discrete complex eigenvalues Xj = cij + ibj, 
Xj = cij — ibj and eigenvectors Wj = uj + ivj, Wj = Uj — ivj, with j = 1, 2, N. 
Thus, iti, t>i, un, vn form a basis of the space 1Z 2N and the invertible matrix 
V = [vi, Mi, , Vn, un] leads us to the Jacobi normal form: 



B = V~ X AV = diag 



aj -bj 



bj aj 



of a 2N x 2N matric B with 2x2 blocks along its diagonal. In this way, using 
the transformation z = V~ 1 x, we can reduce the initial problem to a system 
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Fig. 7. (a) GALIs for 5 coupled standard maps with 
Ki = -2.30, K 2 = -2.35, K 3 = -2.40, # 4 = -2.45, K 5 = -2.50), B = 0.00001 and 
initial conditions that excite 5 different normal modes. Note that the GALI^ for 
k = 2, 3,4, 5 fluctuate around a non-zero value implying a regular motion that lies 
on a 5D torus, (b) The energies for the corresponding normal modes. 



of uncoupled equations, whose evolution is described by the equations: 

z ' = V^AVz, (13) 

where z' = k^ +1 , k% +1 ) T , z = k*, 1%, k^) T and z represents 
the linear normal modes of the map. Thus, in order to excite a continuation 
of one or more of these modes, we extract the appropriate x-initial condition 
by the transformation x = Vz. The evolution in the uncoupled coordinates is 
given by the transformation: 

l"n+l = a ^n ~ bik 3 n (14) 
^n+i = bil J n + a>ik 3 n 

where each pair of (P , /c J ) corresponds to a normal mode of the system. Since 
the map is area-preserving a| + 6| = 1 and we may define the quantity E 3 n = 
(^n) 2 + (^n) 2 as the energy of each mode, which is preserved under the evolution 
of the linear map, i.e. E 3 n+1 = E J n . 

Thus, to study recurrences, we start with the example of iV = 5 coupled 
standard maps, choose a 'small' coupling parameter (3 = 0.00001 and excite 5 
different modes. As a result, we are able to detect a quasiperiodic orbit that 
lies on a 5D torus, indicated by all the GALIfc, k = 1,...,5, approaching a 
constant, as shown in Fig. [71 This is verified in Fig. [TJd, where we plot the 
corresponding energy modes of the coupled system and find that they are 
nearly invariant, as they would have been in the uncoupled case. 

Next, we increase the number of coupled maps from 5 to 20, keeping (3 = 
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Fig. 8. (a) Exciting only one mode in a case of 20 coupled standard maps, with 
(3 = 0.00001 and K j in triplets of -1.35,-1.45,-1.55, the GALI fc , for k = 2,. ..,6 
asymptotically become constant, implying that the motion lies on a 6D torus, (b) 
An orbit with the same initial condition and Kj as in (a) but with a larger cou- 
pling parameter, (f3 = 0.01), becomes chaotic, as all the GALIs are seen to decay 
exponentially to zero. 



0.00001 and choosing different Kj (in triplets of -1.35, -1.45, -1.55), with j = 
1, 20. Exciting now only one normal mode (the 21 st ), much like the original 
FPU case, we calculate the GALIs again and find that the orbit is again 
quasiperiodic, lying on a 6-dimensional torus, as evidenced by the fact that 
GALIfc for k = 2, ...,6 fluctuate around non-zero values, see Fig. [8^l. This is, 
of course, a low-dimensional torus, since generically one expects regular orbits 
of this system to lie on N = 20-dimensional tori! From the calculation of the 
normal mode energies in this experiment, we observe that the initial energy 
is now shared by 6 modes and executes recurrences, for as long as we have 
integrated the equations of motion. The energies of these 6 modes are not 
shown here since their values are too close to each other to be distinguished 
in a graphical representation. 

Finally, as our last experiment, we increase the coupling parameter to j3 — 
0.01, for the same choice of initial conditions and Kj values as in Fig. [8k The 
result is that the influence of the coupling is now quite strong and leads to the 
breakdown of recurrences and the onset of chaotic behavior. This is clearly 
depicted in Fig. [8]d, where the GALI^s of this case exhibit exponential decay, 
implying the chaotic nature of this orbit. 

5 Conclusions 

Localized oscillations in one-dimensional nonlinear Hamiltonian lattices con- 
stitute, for the last 15 years, one of the most active areas of research in Math- 
ematical Physics. Of primary interest has been the discovery of certain exact 
periodic solutions in such lattices, called discrete breathers, which are exponen- 
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tially localized in real space. It is known that when these solutions are stable, 
there are regions around them in phase space, where orbits oscillate quasiperi- 
odically for very long times, even though the presence of linear dispersion is 
expected to lead them to derealization away from the discrete breather. How 
long do they stay, however, in that vicinity and what is the dimensionality of 
the tori on which they lie? Furthermore, how large are these regular regions 
and at what initial conditions (or parameter values) does chaotic behavior 
begin to arise? 

Similar questions can also be posed about another form of localization in 
Fourier space, which has been very recently proposed as an explanation of 
the famous recurrence phenomena known to characterize the dynamics near 
linear normal modes of nonlinear lattices, since the early 1950s. Clearly, the 
interest in such localization lies in the (in) ability of Hamiltonian lattices to 
exhibit energy propagation and equipartition as expected from the theory of 
statistical mechanics. 

To answer such questions, we have applied in this paper a new method that our 
group has recently developed, which is able to distinguish order from chaos, 
identify the dimensionality of tori and predict slow diffusion in thin chaotic 
layers, faster and more reliably than other methods (see also our recent publi- 
cation [13]). This is accomplished by computing what we call the GALI indices 
of the orbits under study, for which we have obtained analytical formulas that 
are in excellent agreement with our numerical calculations. 

We have thus considered two examples of high dimensional systems: One 
Hamiltonian lattice described by 31 Newton's second order differential equa- 
tions and one system of 20 coupled area-preserving pairs of difference equa- 
tions mapping the plane onto itself. In each of these models, we have been able 
to demonstrate the usefulness of the GALI indices in analyzing the properties 
of regular motion localized in real and Fourier space. Furthermore, we have 
shown that our indices can accurately determine cases where the motion be- 
comes delocalized and starts to diffuse chaotically in phase space, long before 
this is observed in the actual oscillations. 

We wish to emphasize the superiority of our approach to the traditional calcu- 
lation of Lyapunov exponents, still used by many researchers to distinguish reg- 
ular from chaotic motion. As we have shown, Lyapunov exponents can be mis- 
leading, as they often tend to decrease in time giving the impression that the 
motion is quasiperiodic, while it is in fact weakly chaotic. The GALI indices, 
on the other hand, are computationally more efficient, as they concentrate on 
the tangent space of the dynamics and explore the linear (in) dependence of k 
unit deviation vectors about an orbit representing the volume of their asso- 
ciated parallelepiped. If, for k = 1,2..., d, that volume remains constant and 
decays for k > d by well-defined power laws, the motion is quasiperiodic on a 
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d-dimensional torus, while if the volume decreases exponentially for all k the 
motion is chaotic. We, therefore, hope that the virtues of the GALI method 
emphasized in this paper, will encourage researchers studying the dynamics of 
multi-dimensional conservative systems to use them in order to study many 
other physically interesting applications. 
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